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FOREWORD 


The Jou I e-Thomson Expander and Heat Exchanger Program was performed 
under NASA Contract NAS2-8984, by the AiResearch Manufacturing Company of 
California, a division of The Garrett Corporation, This work was sponsored 
by the Ames ^Research Center, Moffett Field, California, and administered 
under the direction of Mr. John Vorreiter, M/S 244-14. 

This report was prepared by Mr. R. H. Norman, who was responsible for the 
analytical effort. Mr . J. R, Wenker was responsible for the development and 
fabrication of the assemblies, under the Program direction of Mr. 0. A. Buchmann. 
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SECTION 1 


INTRODUCTION AND SUMMARY 


INTRODUCTION 


The Jo u I e-Thomson Expander and Heat Exchanger Program provides for the 
design, fabrication, and delivery of three assemblies, each consisting of the 
fol lowing components; the combination thereof referred, to as a JTX: 

o I n I et f i I ter 

• Counter flow heat exchanger 

• Jou I e-Thomson expansion device 

• Low-pressure jacket 

The program objective was to develop a JTX, which, when coupled to an open-cycle 
supercritical helium refrigerating system (i.e., a storage vessel), would 
supply superfluld helium (He II) at 2 K or less, for cooling Infrared detectors. 

PROGRAM REQUIREMENTS 

A four-task program was originally planned, including two additional design 
Iterations based on testing by the Government. The first task was design of 
a JTX assembly to provide 30 milliwatts of cooling at 2 K. The second task 
was a performance analysis (prediction) of the JTX over a range of helium 
Inlet temperatures and pressures. The third task was fabrication and delivery 
of an assembly to the Government for testing. The second JTX design was to 
await the test results of the first; similarly, the third design was to await 
the test results of the second. The fourth task was submittal of a final tech- 
n.Ical report after the two design Iterations were completed. 

These requirements were somewhat modified since the Inception of the 
program. All three tasks were carried out for JTX No. 1. JTX No. 2 has been 
fabricated and delivered', but not tested. The design of JTX No. 3 was 
redirected to satisfy the requirements of an astronomical detector cooling 
system, as specified by an amendment to this contract. 

SUMMARY 


All three cryostats were designed; the first, JTX No.. 1, was tested and 
returned for modifications. This modified assembly, designated JTX No. 1A, 
was then delivered to the (Government. The. second cryostat was fabricated, flow- 
tested with nitrogen, and delivered; the third has been designed. 
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The design of all units was based on existing technology developed for 
nitrogen J-T cryostats, with a minor modification to the expansion device. The 
primary design tool for He II cryostats is a JTX Performance Prediction Program. 
Test data has essentially verified the computer model although the capillary 
tube characterization still is difficult for certain input conditions. The 
computer program was modified several times as flow calibration data became 
available, but the basic elements remain the same for JTX No. 3 as they were 
for JTX No. 1 . 

The test data for JTX No. 1 successfully demonstrated the formation of 
superfluid He II at 1.85 K. The test instrumentation, however, did not enable 
a direct comparison of heat exchanger performance to computer predictions except 
for flow rate and cooling capacity. The test flow rate was about 50 percent 
higher than predicted. There are several possible explanations for this: one 

is the tolerance in capillary tube diameter (a 20 percent increase above nominal 
could give 50 percent more flow); another is the possible formation of superfluid 
in the capillary which is not modeled in the computer program. 

JTX No. 2 was to be Identical to JTX No. 1 except for replacing the smal I 
capillary tube with another of the same diameter having twice the length. At 
this point it was established that the original 0.0508 mm (0.002 In.) I.D. 
capillary tubing was no longer available (except by special mill run). A sub- 
stitute of 0.102 mm (0.004 In.) I.D. tubing was approved. The unit was fabri- 
cated, nitrogen flowtested, and delivered. The room-temperature nitrogen 
flow test showed that JTX No. 2 had about half the flow capacity of JTX No. 1, 
which was the desired effect. 

The modified JTX No. lA was to have the same capillary tube replacement 
as JTX No. 2. In the absence of 0.0508 mm I.D. tubing, 0.0838 mm (0.0033 in.) . 
tubing was used. A nitrogen flow test showed JTX No. 1A to have about half 
the flow capacity of JTX No. 2. 

JTX No. 3 was designed for use in the optical cavity dewar of an Astro- 
nomical. Detector Cooling System, for which a preliminary design was conducted 
under Contract Amendment No. 1. This system design was completed and documented 
in a separate report (Reference 1). JTX No. 3 was designed to produce 100 
milliwatts of cooling at 2 K with inlet conditions of 13.6 atm (200 psla) and 
17 K. 
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SECTION 2 


DESIGN AND FABRICATION 


JTX NO. 1 


The JTX assembly consists of the following components: 

• Mandrel andclosure 

• Sheath ( low-pressure jacket) 

• High-pressure tubing, partly finned 

• Expander (2 capillary tubes in series) 

• Fin seals 

• i n 1 et f i i ter 

• GFE f I anges 

The JTX design is derived from existing technology developed for nitrogen 
Jou I e-Thomson cryostats. All of the components listed above are essentially 
the same as a nitrogen J-T' cryostat with the exception of the J-T expander. 

With nitrogen cryostats operating between 150 atm and 1 atm, a single expansion 
capillary is used. To expand supercritical helium (5 to 20 atm) to the low 
pressures required to produce superfiuid helium II (2 to 20 torr, 0.0026 to 
0,026 atm) in a single capillary tube would result in extremely high exit 
velocities and relatively high exit-expansion losses. It is believed that 
a two-step expansion process, where a larger capillary tube is attached to 
the end of a smaller one, promotes a more controlled expansion into the vacuum 
of the cooling cavity. 

The JTX assembly, without flanges and filter, is presented in Figure 
2-1, showing the approximate dimensions for JTX No. 1. The mandrel is a 
seamless tube of L605 (Haynes Alloy 25), a high cobalt alloy with desirably 
low thermal conductivity and heat capacity. The mandrel Is sealed at the 
cold end to prevent helium flow from bypassing the finned tube heat exchanger. 
The sheath is also fabricated of L605 alloy thin-wall tubing. The diameter 
is selected for proper mandrel— to— sheath gap for the selected fin. The copper 
fin is formed from 0.203 mm (0.008 in.) diameter wire to dimensions of approxi- 
mately 0.104 rnn (0.0041 in.) by 0.318 mm (0.0125 in.) and annealed. A post 
heat-treat cleaning is performed prior to winding the fin. The high-pressure 
tubing of 300-series stainless steel of 0.559 mm' (0.022 in.) O.D. by 0.305 
mm (0.012 in.) I.D. is cut to extra length and cleaned to remove contaminants 
which could react in the brazing cycle. Following internal cleaning the tubing 
is crimped shut and externally cleaned and plated. Copper and silver are plated 
onto the high-pressure tube to provide braze material for the silver-copper 
eutectic which bonds the fin to the tube. The copper fin is then edge-wound 
over the high pressure tube, with a slight excess temporarily tack-soldered 
at each end to hold it in place. The fin spacing is control led to about 39 
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1.2 IN. 

Figure 2-1. Design Details of JTX No. 1 (Without Flanges) 


fins/cm (100 fins/in.). The wound finned tube is Then brazed in an inerT 
atmosphere to secure the fin to the tube with a thermally conductive joint. 

A careful visual inspection of the braze fillets is made to assure that a 
proper f in-to-tube bond exists. A short length of unfinned tube is wound 
around the mandrel and tack-brazed at each end. The high-pressure tube is 
then cut to reopen it, and the first section of capillary tube is inserted 
a short distance in and sealed with a braze joint. This capillary tube is of 
300-series stainless steel with nominal dimensions of 0.153 mm (0.006 in.) 

O.D. by 0.0508 mm (0.002 in.) I .0.* by 7.4 cm (2.9 in.) long. The second 
capillary tube is slipped over the first, sealed with a braze, and the two 
are wrapped onto the mandrel and secured with soft solder. The nominal dimen- 
sions of this second capillary are 0.330 mm (0.013 In.) O.D. by 0.178 mm (0.007 
in.) l.D. by 22 cm (8.5 in.) long. The finned tube and expansion capillaries 
are then wrapped around the mandrel. The inner fin seal, 0.356 mm (0.014 in.) 
diameter nylon, is installed and adjusted to provide the proper fin pitch 
and outside diameter. The nominal O.D. of the fin is 1.25 mm (0.049 in.), 
therefore approximately 20 coils provides the required 40 cm (15.8 In.) length 
of finned tube, occupying about 2.5 cm (1.0 in.) of mandrel length. The outer 
fin seal, 0.406 mm (0.016 in.) diameter nylon, is then installed and secured. 

The wrapped heat exchanger section is then carefully sized using tooling to 
ensure that the wrap is uniform and will properly contact the sheath. The 
fin seals are used to direct the low pressure helium away from the sheath 
and mandrel and into better contact with the finned tube. 

The photograph in Figure 2-2 shows the final assembly of heat exchanger 
and expansion tube wound around the mandrel. The sheath with GFE flanges 
and inlet heater are shown beside the heat exchanger/expander subassembly in 
Figure 2-3. The large upper flange contains the instrumentation feedthrough 
provisions and the lower smaller flange mates to the sensor mounting plate. 

The helically coiled tube between the flanges contains the JTX inlet electrical 
heater. The photograph of Figure 2-4 shows the completed JTX assembly, prior 
to installation of sensors and instrumentation. The soldered sheath penetration 
of the inlet tube can be seen just below the helical heater. Above the large 
flange is the^ inlet filter assembly fitted to a 0.318 cm (1/8 in.) O.D. supply 
line. This filter is of the stainless steel sintered-mesh type with a nominal 
2 micron rating. 

The modified JTX assembly, JTX No. 1A, was identical in design to JTX No. 1 
with the exception of the expansion capillary tubes. A small capillary tube 
which measured (single microsection) 0.185 mm (0.0073 in.) O.D. by 0.0838 
mm (0.0033 in.) l.D. by 60 cm (23.6 In.) long, was installed. The larger 
capillary tube was replaced by a 30 cm (12 in.) length of 0.457 mm (0.018 in.) 
O.D. by 0.254 mm (0.010 in.) l.D. tubing (not measured). 

JTX NO. 2 

The second assembly, JTX No. 2, was identical to JTX No. 1, except for 
the expander. In this unit a small capillary of nominal dimensions 0.203 mm 


^An approximate measurement of this tubing from a 400X photograph revealed 0.007 
in. O.D. by 0.0024 to 0.0026 in. l.D. 
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Figure 2-2. J-T Expander and Heat Exchanger Wound on Mandrel 
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76403-2 


Figure 2-3. J-T Expander and Heat Exchanger with 
Sheath, Flanges and Inlet Heater 
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Figure 2-4. Completed JTX Assembly 
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(0.008 In.) O.D. by 0.102 mm (0.004 In.) I.D. by 86 cm (34 In.) long was 
Installed. The large capillary tube was the same size as that used on JTX 
No. 1A. The actual measurements from a microsection of the small capillary 
tube were 0.0081-0.0082 in. O.D by 0.0044-0.0046 In. I.D. 

The sheath, flanges, and Inlet filter were assembled around the wrapped 
mandrel, as with JTX No. 1, but the helically colled heater was omitted. An 
extended length of inlet tubing was provided between the upper flange and 
the point of sheath penetration. 

JTX NO. 3 


The third JTX assembly was designed to provide cooling for Infrared 
detectors in an optical cavity dewar of the Astronomical Detector Cooling 
System (Reference 1). The design Is essentially the same as that for the 
preceding two cryostats with the following exceptions. The finned tubing 
of the heat exchanger is longer: the total length Is 61 cm (24 In.), occupying 

about 4.1 cm (1.6 In.) of mandrel length. The expander capillary tubes are 
the same nominal sizes that were used on JTX No. 2, except that the length 
of the smaller tube Is 46 cm (18 in.). The flanges, fittings. Inlet filter, 
and sheath assembly are all different than in previous JTX's in order to be 
Installed in the relatively small envelope of the optical cavity dewar. These 
details are shown in a layout drawing presented In Reference 1. 

A summary of expander capillary tube sizes for all four assemblies is 
presented in Table 2-1. 

FLOW CALIBRATION TESTS 

After final assembly and cleaning of each JTX assembly, a flow test with 
ambient (room) temperature nitrogen and/or helium is performed. High-pressure 
gas is regulated Into the JTX and vented through the expander, exhausting to 
atmospheric pressure. The flow rate is measured as a function of inlet 
pressure. At the moderate pressures used, the nitrogen cools very little upon 
J-T expansion. The helium, of course, warms Insignificantly upon J-T 
expansion at room temperature. 

These preliminary tests are performed to verify that the expander capillary 
is not plugged. In the production of nitrogen J-T cryostats for other applica- 
tions this test Is routinely performed to "calibrate” the capillary tube; 
that is, the required flow range is obtained by clipping off the capillary 
tube and re-testlng. In the JTX application this calibration cannot be performed 
because of the lack of a cryogenic supply of supercritical helium and a special 
test setup. If liquid nitrogen were used, it would not verify He II performance 
without actually prior -testing with supercritical helium, anyway. 

The JTX No. 1 assembly was flow-tested with helium and nitrogen. The 
other two assemblies, JTX No. 2 and JTX No. 1A were tested with nitrogen 
only. The four sets of data points are plotted in Figure 2-5, with the test 
data correlations obtained by linear regression analysis (a power curve fit 
was chosen, since the data plotted on a log-log scale indicated a straight 
line). 
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TABLE 2-1 


SUMMARY OF EXPANDER CAPILLARY TUBE SIZES 


- 

JTX Assembly 

Dimension 

No. 1 

No. lA 

No. 2 

No. 3* 

Sma i 1 Capi 1 lary Tube 





Nominal I.D., mm 

0.0508 

? 

0.102 

0.102 

in. 

0.002 


0.004 

0.004 

Nominal O.D., mm 

0.153 

? 

0.203 

0.203 

In. 

0.006 


0.008 

0.008 

Measured I.D,, mm 

0.0635 

0.0838 

0.114 



In. 

0.0025** 

0.0033 

0.0045 

— 

Measured O.D., mm 

0.178 . 

0.185 

0.208 



in. 

0.0070** 

0.0073 

0.0082 

— 

Length, cm 

7.4 

60 

86 

46 

in. 

2.9 

23.6 

34 

18 

Large Capi 1 lary Tube 

■1 




Nominal i.O., mm 


0.254 

0.254 

0.254 

in. 


0.010 

0.010 

0.010 

Nominal O.D., mm 

0.330 

0.457 

0.457 

0.457 

in. 

.0.013 

0.018 

0.018 

0.018 

Length, cm 

22 

30 

30 

30 

in. 

8.5 

12 

12 

12 


*These sizes are design values used 


in performance prediction analysis 


^^Approximate measurement from photo of 


contaminated section 
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Figure 2~5. Correlation of Flow Calibration Test Data 



In addition to the expander flow tests, a special nitrogen test was 
performed with JTX No. 2 to measure the pressure drop through the' low-pressure 
side finned heat exchanger as a function of flow rate. This was accomplished 
by blocking the sheath on the expander end and attaching an inclined manometer 
to measure the pressure differential between the expander outlet and the atmos- 
phere (the flow being exhausted out the sheath/mahdre I annulus at the JTX 
inlet end). The purpose of this special test was to correlate the Zeta-factor 
in the analytical model to actual pressure-drop data. Again, the validity 
is limited by the- measurements being taken with nitrogen at room temperature 
and atmospheric pressure, rather than helium at 4 to 20 K and subatmospher Ic 
pressure. An important result of this test was, nevertheless, the modification 
of the computer model to account for 'Maminar” flow around the fins, for very 
low Reynolds numbers. This will be discussed in the next section. Figure 
2-6 shows the good correlation between the pressure-drop data and the laminar- 
flow-modified computer model, for a Zeta-factor of 0.748. This Zeta was obtained 
by reducing the pressure-drop data to Zeta-f actors, and computing the arithmetic 
mean value (the standard deviation was +0.010). 
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PERFORMANCE ANALYSIS 



SECTION 3 


PERFORMANCE ANALYSIS 


J-T THERMODYNAMICS 


The Jou I e-Thomson (J-T) process can best be understood by examining a 
pressure-enthalpy diagram of the working fluid. Figure 3-1 presents a schematic 
of the J-T expander and heat exchanger, together with a gross representation 
of the pressure-enthalpy diagram for helium 4 In the low-temperature region. 

The diagram shows the state points corresponding to the key physical locations 
In the JTX as depicted In the. schematic. 

Briefly, the process consists of two heat exchange steps In which high- 
pressure gas Is cooled (1 to 2) with low-pressure vapor (4 to 5) boiling out 
of the two-phase mixture (3); separated by an Isenthalpic (adiabatic Irreversible] 
expansion to a low pressure where coexisting liquid (He II, in this case) and 
vapor (He I) appear. The refrigeration load (IR sensor) and heat leak cause 
the vapor to boll off, returning to the heat exchanger, leaving behind the liquid 
phase to maintain a stable temperature level. The thermodynamics of this process 
imply that the formation of superfluid helium (He II) Is possible if suitable 
expander inlet temperatures and chamber back pressures (vacuum level) are obtalnec 

COMPUTER PROGRAM DEVELOPMENT 


A computer program, written In FORTRAN V language for the Uni vac 1108 
computer (University Computing Company software). Is used to predict the per- 
formance of a JTX of a given design (all geometrical parameters must be specified] 
The only operational parameters required are the heat exchanger hlgh-pressure-sld< 
inlet pressure and temperature, and the refrigeration load temperature (which 
determines the helium pressure at the Inlet to the low-pressure side of the 
counterflow heat exchanger). Helium thermophysical properties maps were assemblet 
to cover the expected operating pressure and temperature ranges. The overall 
logic and operation of the computer program are discussed In this section, 
while a detailed summary of the key subroutines Is presented In Appendix A. The 
development of a pressure-enthalpy diagram extending into the helium II region, 
and a tabulation all properties required by the program are presented In Appendix 

A block diagram of the JTX computer program showing the subroutine structure 
Is presented in Figure 3-2. There are four basic groups of routines under the 
main program (MAIN): (1) a data element containing the helium thermophysical 

properties (HELIUM/DATA) — specific enthalpy, specific volume, dynamic viscosity 
and thermal conductivity, each as a function of pressure and temperature, (2) 
a group of service routines for reading and looking up the helium properties 
(HEREAD, THERMO, etc.), (3) the main operational routines which perform the 
J-T cycle analysis (JTX) and characterize the pressure-drop and heat transfer 
performance (DPIHX, EXPDR, HEATX, etc), and (4) the purely mathematical Iteration 
routine NEWTON, and Its own service subroutines (LNMOD, MATMUL, etc), which 
performs the iteration on the two key system variables from which all other 
unknowns are directly calculated. The subroutine FENCE contains the limits 
imposed on the two independent variables converged upon by NEWTON. 
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The two key variables (unknown) are expander inlet temperature and flow 
rate. The MAIN program provides a guess of each one and then calls NEWTON. 

This proprietary AiResearch program makes use of the so-called Newton-Jacobian 
matrix method of solving an array of unknowns by imposing stepwise incremental 
changes and examining the derivative of each function of that variable to 
obtain the next guess in the path toward convergence. NEWTON calls the sub- 
routine JTX every time it obtains a new guess. 

Subroutine JTX is the cycle analysis program which calculates the system 
.pressure drops and heat exchanger performance, based on the guessed value of 
expander inlet temperature and flow rate. It starts from the high-pressure-side 
inlet conditions and cails DPIHX to compute pressure drop inside the coiled 
tube. This pressure drop subtracted from the JTX inlet pressure gives the 
Inlet pressure to the expander. With the expander Inlet conditions (P, T) 
known, EXPDR is called to calculate the pressure drop in the two capillary tubes, 
his is the most difficult calculation in the program because of the widely 
varying fluid properties from inlet to outlet. The procedure is to divide 
each tube into 20 increments of length and calculate the isenthalpic pressure 
drop in each segment by the conventional incompressible flow equation, assuming 
fluid properties are constant over the length increment. The friction factor ' 
correlation used in both DPIHX and EXPDR is taken from Reference 2, for turbulent 
flow in coi led tubing. 

The pressure drop on the low-pressure side of the heat exchanger Is 
computed by DPOHX, from the inlet conditions given by the cooling cavity. 
Saturated vapor at the refrigeration temperature enters the annulus and flows 
oyer the finned tubing. The friction factor correlation for flow around finned 
circular tubes is based on Reference 3, with the addition of a laminar-type 
flow dependence at very low Reynolds numbers. This subroutine makes use of 
the Zeta-factor in the calculation of flow area. This factor is determined 
from experimental data (see Figure 2-6), and attempts to account for the flow 
restriction caused by the presence of the fin seals, as well as the effect 
of continual expansion and contraction as the helium flows around the tube 
coils, fin seals, and between the fins. 


The heat exchanger performance is calculated by subroutine HEATX, given 
all the inlet and outlet conditions. At this point, the outlet temperature on 
the low-pressure side is obtained by heat balance over the heat exchanger. 

The output of HEATX is the overall thermal conductance, based on a coiled-tube 
heat transfer coefficient (inside tube) correlation from Reference 4, and a 
correlation for the outside finned tube coefficient from Reference 3. The 
overall conductance (UA) is then used to calculate the heat exchanger effective- 
ness, which in turn gives the low-pressure-side outlet temperature. This 
temperature, T5 in Figure 3-1, is one of the convergence parameters. When the 
I5 obtained from HEATX agrees with the T5 obtained from the heat balance, then 
one of the convergence closures is satisfied. The other convergence closure 
is the expander pressure drop (APx), which is made to agree with the difference 
between the state point pressures P2 and P3 in Figure 3-1. In the development 
of the computer program, this latter convergence closure became the focus of 
much attention. In many early cases this closure failed to converge, while the 
T5 convergence was excellent, and the cycle state points all seemed reasonable. 
Detailed study of the subroutine EXPDR (a stand-alone version was developed) 
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revealed that the expander pressure drop was extremely sensitive to flow rate, 
and the closer the flow rate was to the converged value, the more sensitive it 
became. Convergence was finally obtained for most cases by comparing the logarithm 
of the pressure differences, instead of the pressure differences themselves. 

The resulting flow rates were usually within a few percent of the unconverged 
value. This problem is believed to be caused by the iarge change in fluid 
properties (density and viscosity) over the range of pressures encountered in 
the expander; from supercritical, high-density fluid, passing very near the 
critical point, then into the two phase region where liquid and vapor coexist. 

At first, a conservative approach was taken by assuming that the viscosity 
of the two-phase mixture was the viscosity of the vapor only. This apparently 
created large discontinuities between two successive length increments, making 
convergence more difficult. To obtain a higher percentage of converged results, 
a mixed-mean or "homogenized”, viscosity value was used. The mixed value is 
simply the mass-weighted average of the liquid and vapor values. 

As will be shown later, this model predicts JTX performance reasonably 
well. The same logic was used in another computer program for predicting the 
performance of nitrogen J-T cryostats. The predicted flow rate agreed very 
well with measured flow rates at the normal operating points. These nitrogen 
cryostats were part of a production order where each unit has its expander 
capillary tube "calibrated" by clipping off and running a flow test. At first, 
predictions of required capillary tube length were low by a factor of two. 

After measuring the inside diameter of a sample of the capillary tubing 
(same tubing as used for the small capillary in JTX No. 2), the computer pre- 
dictions based on the measured diameter fell within the measured flow rates 
at the actual tube lengths used. Figure 3-3 presents the results of this 
capi I I ary tube study. 

The service subroutines which read and perform Interpolation of the helium 
thermophysical properties are all quite specialized. The helium data are re- 
stricted to the pressure and temperature ranges within which the JTX is 
typically expected to operate. The lookup options are controlled by the 
particular subroutine calling the lookup subroutine. Table 3-1 presents a brief 
description of the helium properties subroutines and their limitations. A com- 
plete tabulation of the thermophysicai properties used in the program is given 
in Appendix B. 

A sample printout of the JTX performance prediction program is presented 
in Figure 3-4. The format is identical for each case, regardless of the number 
of input changes made. All input parameters are printed out for each case. 

"Design Parameters" gives the geometrical Input, and heat leak, if known; and 
"Input Conditions” gives the three state points required to define the operat- 
ing conditions of the case. Every attempt was made to give self-descriptive 
parameter names. The cycle state point identification numbers refer to those 
given in Figure 3-1. 
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Figure 3-3. Nitrogen J-T Cryostat Capillary Tube Study 
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TABLE 3-1 


DESCRIPTION OF HELIUM PROPERTIES ROUTINES 


Name 

Description 

HELIUM/DATA 

Data element containing helium single-phase properties 
as a function of temperature and pressure. The tempera- 
ture range Is 2.222-20 K (4.0-36 R). The pressure ranges 
are as fol lows; 


Specific enthalpy 0-20.41 atm (0-300 psia) 


Specific volume and 0.068-20.41 atm 

dynamic viscosity (1-300 psia) 


Thermal conductivity 2.244-20.41 atm (32.99-300 psia) 

REREAD 

Reads In helium single-phase thermophysical properties 
from data element HELIUM/DATA, and converts data arrays 
to MAPRDY format. 

SATPRP 

Helium saturation properties lookup by linear Interpola- 
tion. Contains 41 -point data for specific enthalpy, 
specific volume, and dynamic viscosity from 0.833 K 
(1.5 R), taken as zero pressure, to the critical point 
5.202 K (9.363 R), including the lambda point 2.177 K 
(3.919 R). Saturated liquid and vapor properties and 
saturation pressure are given at each temperature. 
Saturated liquid viscosity below the lambda point 
(He 11} is given as Infinitesimally small number since 
It is not used In program. Four lookup options;- 


(1) Pressure in, sat. liquid props, out 
I (2) Pressure In, sat. vapor props, out 

(3) Enthalpy In, sat. liquid props, out 

(4) Temperature In, sat. liquid props, out 

THERMO 

Three-dimensional lookup of helium single-phase 
specific enthalpy and specific volume .read Into MAIN 
program by subroutine HEREAD. Uses MAPRDY for 3-D 
lookups. Calls SATPRP for saturation data when input 
P and T fall between single-phase data point and two- 
phase dome. Linear Intepolation then employed. Two 
lookup options; 

- 

(1) Pressure and temperature In, enthalpy and 
specific volume out 


(2) Pressure and enthalpy In, temperature and 
specific volume out 
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TABLE 3-1 (Cont) 


Name 

Description 

THERMO 

(Cont) 

With option (2), if pressure is less than critical 
and enthalpy puts lookup point just outside dome, 
then it must be on liquid side only. (This only 
occurs in EXPDR, which must have inlet enthalpy less 
than critical or He II may not be produced). Low- 
pressure specific volume (below 1.0 psia) Is 
calculated from ideal gas equation. Specific 
volume data below critical pressure in single-phase 
region is limited to liquid side, also. 

TRANSP 

Three-dimensional lookup of helium single-phase 
viscosity and thermal conductivity read into MAIN 
program by subroutine REREAD. Uses MAPRDY for 3-D 
lookups. Viscosity P-T lookups below critical 
pressure performed as in THERMO, being limited to 
liquid side. Thermal conductivity P-t lookups 
limited to supercritical pressure (HEATX is the 
only subroutine using thermal conductivity, and 
J-T heat exchanger does not operate In inter- 
mediate pressure region). Ail other pressures 
(below critical) use temperature-dependent built-in 
data and linear interpolation, since low-pressure 
transport properties are nearly independent of 
pressure. 

LAGIN2 

Performs linear interpolation and extrapolation of 
any tabular data, using the binary search method of 
reaching the two data points between which the answer 
lies (very fast). 

MAPROY 

Three-dimensional lookup routine used by THERMO and 
TRANSP to find properties as a function of two input 
variables. 

KFIN 

Data and lookup routine (linear interpolation) for 
stainless steel tubing and OFHC copper fins. Built-in 
data provided for temperatures up to 50 K (90 R). 
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Figure 3-4. Sample Computer Printout of JTX Performance Prediction 
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JTX NO. 1 PERFORMANCE PREDICTION 

The computer program was used to predict the performance of JTX No. 1 
oyer a wide range of inlet temperature and pressure, for the design parameters 
given in Figure 3-4, with one exception. At the time these predictions were 
made, the inside diameter of expansion capillary tube No. 1 had not been meas- 
ured, thus the nominal value of 0.0508 mm {0.002 in.) was used. Five inlet 
pressures ranging from 3.4 atm (50 psia) to 20.4 atm (300 psia) were studied 
for a refrigeration load temperature of 2.0 K. After more than 50 runs were 
made, the range of inlet temperature over which the JTX would produce liquid 
II bscame quite clear. The temperature limits of operation are best illus- 
trated on a plot of expander inlet temperature (T 2 in Figure 3-1) versus JTX 
inlet temperature (T| in Figure 3-1), as shown in Figure' 3-5. The lower limit 
of JTX inlet temperature is not a physical limit, rather a helium properties 
limit in the computer program. The program has no enthalpy, density, or 
viscosity data at high pressure (above saturation) below a temperature of 
2.22 K (just above lower lambda point), thus preventing characterization of 
the expander pressure drop and flow rate. No. data for' He II in this region 
could be found. ^ The upper limit of expander inlet temperature is based on a 
value of specific enthalpy Just below the critical-point value. Examination 
of the pressure-enthalpy diagram presented in Appendix B shows that a slightly 
higher enthalpy value may allow He II production at 2 K, although lower tem- 
perature refr igeration would be limited. Even at the critical-point enthalpy 
the amount of liquid produced at 2 K is very small (about 16 percent, by 
weight). Other factors considered in this choice of expander specific 
enthalpy limit are: (1) the uncertainty in the pressure-enthalpy diagram 
for He II, and (2) the possibility of difficulty in obtaining continuous 
values of fluid properties when interpolating on data near the critical 
point (this difficulty is probably the major cause of not obtaining con- 
vergence in the expander pressure-drop/f low rate iteration with the present 
enthalpy I imit) . 

The upper limit of JTX inlet temperature cannot be determined precisely 
because a heat exchanger effectiveness of exactly 100 percent will not con- 
verge in the program. The curves for 3.4 atm (50 psia) and 6.8 atm (100 psia) 
are stopped at inlet temperatures which gave heat exchanger efficiencies of 
99.8 and 99.2 percent, respectively. At the other three pressures the heat 
exchanger efficiencies were 98.3 percent, or less, at the JTX inlet tempera- 
ture of 20 K. The JTX will no doubt operate at temperatures above this, 
although at greatly reduced refr igerat ion capacity. The 20 K limit is the 
upper limit of helium properties built into the computer program. 

Figure 3-6 shows the predicted gross refr igerat ion at 2.0 K as a func- 
tion of JTX inlet pressure and temperature. Note the reduction of cooling 
capacity as the inlet temperature is increased for a given inlet pressure. 

At a given Inlet temperature, however, there appears to be an ''optimum” 
pressure above and below which the refrigeration capacity drops off. The 
curves in Figure 3-6 should be considered as approximate, since "smoothing” 
of many isotherms was required to fill out the plot, especially at the lower 
temperature and pressures, where very few computer cases actually converged. 
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Figure 3-5. Prediction of JTX No. 1 Expander Inlet Temperature 
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Figure 3-6. Prediction of JTX No. 1 Refrigeration Performance 
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The best predictions were made in the region above 6.8 atm (100 psia) and 
9-10 K. The dashed curve at the top of Figure 3-6 shows the approximate limit 
set by a minimum allowable expander inlet temperature of 2.22 K, due to lack 
of high-pressure He II data. 

The predicted flow rates for JTX No. 1 as a function of inlet temperature 
and pressure are presented in Figure 3-7; again, the nominal inside diameter 
of the small capillary tube was used. The same computer runs produced the 
data for all three plots (Figures 3-5, 3-6, 3-7). Note that in Figure 3-7 
as JTX inlet pressure Is increased, the flow rate becomes nearly Independent 
of Inlet temperature. This is probably due to the effect of more uniform 
thermophysical properties as more of the expander operates away from the 
critical point. 

JTX NO. 1 TEST DATA ANALYSIS 

in December 1975, the Jet Propulsion Laboratory (JPL) performed tests on 
JTX No. 1, wherein the assembly shown in Figure 2-4 (with a vacuum container 
housing the portion below the Instrumentation flange) was immersed in a liquid 
helium dewar. The inlet heater allowed the warmup of a pressurized (liquid) 
helium stream to. the desired JTX inlet temperature. 

The preliminary test results (steady-state data) are presented in Table 3-2, 
as reported by Dr. L. C. Yang. 


TABLE 3-2 

PRELIMINARY TEST RESULTS OF JTX NO. 1 BY JPL 


JTX Inlet 
Pressure 
(psia) 

JTX Inlet 
Temperature 
(K) 

Sensor Plate 
Temperature 
(K) 

Warm-End 

Sheath 

Temperature 

(K) 

Mass Flow 
Rate 
( Ib/hr) 

Gross 

Refrigeration 

(mW) 

100 

4.4 

1.81 

5.0 

0.0544 

176.8 

100 

6.6 

1.81 

6.2 

0.0488 

150.2 

100 

10.1 

1.85 

8.3 

0.0528 

95.4 

• 100 

12.1 

1.84 

8.6 

0.0532 

51.5 . 

100 

15.2 

1.84 

11 .0 

0.0520 

25.2 

200 

4.5 

1.9 

5.4 

0.0808 

272.6 


It was generally agreed that the demonstration of achieving refrigeration 
at superfluid He II temperatures, in a stable manner, with a J-t cryostat was a 
great success. However, evaluation of the counter flow heat exchanger efficiency 
was impossible because the warm-end sheath temperature apparently did not mea- 
sure the low-pressure-side outlet gas temperature (T5 In the schematic of 
Figure 3-1) well enough. If the values presented in Table 3-2 were used to 
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represent T5 In the thermodynamic cycle analysis, no heat balances around the 
system are obtained, and for the first and last cases a heat exchanger effec- 
tiveness of over 100 percent is' indicated. A simplified analysis of the heat 
transfer from the outlet gas to the sheath near the warm-end sheath temperature 
measurement, for the third case in Table 3-2, showed that the flowing gas could 
be 1.7 K warmer than the sheath at that point (8.3 K). This difference is 
enough to invalidate the heat balance, and prevents a meaningful comparison 
with the predicted heat exchanger effectiveness (that is 98 percent, versus 
78 percent from the 8.3 K data point). 

The refrigeration and flow rate data, however, can be made to compare 
reasonably well with computer program predictions. Figure 3-8 presents a plot 
of the 100-psia test data, showing gross refrigeration and mass flow rate 
plotted as a function of JTX Inlet temperature. The dashed lines are the pre- 
dictions based on the nominal inside diameter of the small capillary tube 
(0.002 in.). The dashed flow rate curve is taken directly from Figure 3-7. 

The correspond i ng refrigeration curve is taken from the computer results which 
were used to construct the performance map of Figure 3-6, hence do not agree 
exactly with the "smoothed" curves, especially at the higher inlet temperatures. 
The computer results all used 2.0 K as the refr igerat ion temperature, rather 
than the 1.81-1.85 K obtained during testing. The effect. of this difference 
Is expected to be very small, however, as the thermophysical properties 
of helium do not vary significantly between these temperature levels (the 
only significant effect would be in the heat exchanger low-pressure-side 
pressure drop, which was not independently evaluated). 

Figure 3-8 shows that the measurements of both refrigeration and flow 
rate were higher than the predicted curves using the nominal capillary tube 
diameter. As noted in Section 2 of this report, the small capillary tube size 
was later measured to be 20 to 30 percent larger in diameter (at one point, 
anyway) than Its nominal value. Therefore, new computer predictions were 
made, using a 25 percent larger inside diameter, and a comparison with the 
test data showed much better agreement. This is illustrated in Figure'3-8 
by the sol id- line curves, which appear to be within experimental uncertainty 
of the test data. The lower inlet-temperature test cases cannot be analyzed 
by the computer model because of the low-temperature, high-pressure. He II 
thermophysical properties limitations discussed previously. The lowest 
temperature tested is lower than the JTX is likely to ever encounter from 
an actual supercritical helium storage tank. Thus, the low-temperature 
expander Inlet temperature limitation of the computer program is not con- 
sidered to be a very serious shortcoming. 

A performance prediction of JTX No. 1A (same as JTX No. 1 except for a 
.longer, larger diameter, small capillary tube) was made for the third case con- . 
ditions of Table 3-2, namely, the case of 10.1 K inlet temperature at 100 psia. 
The results for a refr igeration temperature of 2,0 K were; mass flow rate, 

0.0177 kg/hr (0.039 Ib/hr); and gross refrigeration, 55 milliwatts. A glance 
at the nitrogen flow calibration test data of Figure 2-5 reveals a large 
difference in the predicted flow rate at He II levels as compared with room- 
temperature nitrogen data. The ratio of JTX No, 1 flow to JTX No. 1A flow 
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at the 2.0 K level is predicted to be 1.50, while the same ratio from the 
nrtrogen room-temperature data !s about 4*4. This apparent conflict remains 
unexplained at this time. No other data exists to substantiate the indica- 
tion that room-temperature flow data cannot be correlated to He II flow data, 
as was expected. 

JTX NO. 2 PERFORMANCE PREDICTION 


The preHminary analysis of JTX No. 2 performance was based on the 
original design intent — that is, the same as JTX No. 1 except for doubi ing 
the length of the small capillary tube. This was done, over a narrower range 
of input parameters (27 cases) than for JTX No. 1, before it was established 
that the small diameter tubing was no longer available. The next larger size 
tubing was incorporated after estimating how much coiled length of it'would 
fit on the existing mandrel. This length was estimated to be 86 cm (34 in.), 
assuming a double wrap of 25 coils, followed by 17 colls (12 in.) of the larger 
(0.018 in. O.D.) capillary tube. The performance prediction computer program 
was run for these dimensions and the same inlet conditions as several of the 
JPL test cases.. Table 3-3 presents the flow rate results, showing also the 
predicted values for the capillary tube size originally intended for JTX No. 2 
(with a 20 percent increased inside diameter over the nominal size). Although 
many of these cases were unconverged, the predicted flow rate is probably 
with 5 to 10 percent of the converged value. The conclusion reached here 
was to expect the JTX No. 2 flow rate to be approximately the same as the 
JTX No. 1 test data. Again, as noted above for the case of JTX No. 1A, the 
He II predicted flow performance is at odds with the room- temperature nitrogen 
flow calibration test data presented In Figure 2-5. 

TABLE 3-3 

JTX NO. 2 PREDICTED FLOW RATE 

FOR 1.84 K REFRIGERATION (100 PS I A CASES) 


JTX Inlet 

Predicted Flow Rate (Ib/hr) 

JPL Test 

Temperature 






0.0024 in. 1 .D. 

0.0024 in. !.D. 

0.004 in. I.D. 

r 1 kJW r\Cj 1 e 


2.9 in. long 

5.8 in. long 

34 in. 1 ong 


(K) 

cap. tube 

cap. tube 

cap. tube 

( Ib/hr) 

10.1 

0.0530 

0.0388* 

0.0552* 

0.0528 

12.1 

0.0511 

0.0335 

0.0520* 

0.0532 

15.2 

0.0490* 

0.0331 

0.0485* 

0.0520 


*Unconverged results, due to expander pressure— drop/f low rate problem; 
flow rate probably within 5-10 percent of correct value. 
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JTX NO. 3 PERFORMANCE PREDICTION 


The performance prediction of JTX No. 3 was carried out in the same mannei 
as for JTX No. 1. A total of 48 computer cases were run covering 5 values of 
inlet pressure from 6.8 atm (100 psia) to 20.4 atm (300 psia) and JTX inlet 
temperatures up to 18 K. The selection of the design points and the range of 
performance parameters studied is discussed in Reference 1. 

All results of the computer runs (except for cases outside the minimum 
expander inlet temperature limit) were first plotted as gross refrigeration 
against JTX inlet temperature, for each of the 5 pressures. Unconverged 
results were plotted along with the converged values, and taken into considera- 
tion when drawing the curves, presented in Figure 3-9, Fortunately, the 
percentage of unconverged results was much smal ler over the range of pressures 
and temperatures studied for JTX No. 3. The expander subroutine seems to work 
better at higher flow rates. 

To obtain a better picture of JTX performance as a function of inlet 
pressure, the curves of Figure 3-9 are cross-plotted to obtain the refrigera- 
tion performance map in Figure 3-10. The corresponding predicted flow rates, 
as a function of inlet pressure and temperature, are plotted in Figure 5-11, 
showing the relatively weak dependence on temperature and strong dependence 
on inlet pressure. Both of these figures are reproduced from Reference 1. 
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JTX INLET TEMPERATURE, KELVIN 

Figure 3-9., Predicted JTX No. 3 Refr igerat ion as a Function 
of Iniet Temperature 
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Figure 3-10. .Prediction of JTX No. 3 Refrigeration Performance 
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APPENDIX A 
COMPUTER PROGRAM 



NOMENCLATURE 


Fl ow cross-secti onal area 
A^g Fin surface area 

A. Surface area inside high-pressure tube 

A^ Finned tube surface area on low-pressure side 

Ag Total surface area on low-pressure side 

C Capacity rata ratio = 


'MAX 

Maximum of 

rn 

and Cp^ 


MiN 

Mi ni mum of 

PH 

and Cp^ 


^PH 

Heat capaci ty , 

high-pressure side 


^PL 

Heat capaci ty , 

low-pressure side 


D 

c 

Finned tube coi 

1 diameter 


^ec 

Expander tube 

coi 1 di ameter 



Inside di amete r 

of small expander 

tube 

^e2 

Inside diameter 

of large expander 

tube 


Finned tube overall diameter 
Hydraulic diameter 

Dj Hi gh-pressure tube inside diameter 

Mandrel outside diameter 
High-pressure tube outside diameter 
Sheath inside diameter 

DPV ‘ Volume-normalized low-pressure side pressure drop 
Friction factor (Darcy) in expander tube increment 
j Friction factor (Darcy) inside high-pressure tube 

Friction factor (Fanning) outside finned tube 
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G 


•"c 

H 

h. 

I 

h 

0 

k 


“el 

“e2 


ML^ 

^tu 

P 


ex 

Pr 


AP. 

I 

AP 

o 

Q 

Re 


T 

UA 

V 

W 


Hass velocity = W/A^ 

Constant = 4.1698 x 10® 

Specific enthalpy 

Heat transfer coefficient inside high-pressure tube 
Heat transfer coefficient outside finned tube 
Thermal conductivity 

Expansion loss factor between expander tubes 

Length of small expander tube 

Length of large expander tube 

Total length of finned tube 

Increment length of expander tube 

Fin efficiency parameter 

Number of transfer units = UA/(W* C,,. ,.) 

MIN 

Absolute pressure 
Fi n pi tch 

Expan“der exit pressure 

Prandtl number = Cq p/k 

Pressure drop inside high-pressure tube 

Pressure drop outside finned tube 

Expander pressure drop 

Refrigeration capacity or heat leak 

Reynolds number = D. G/m^ 

h 

Absolute temperature 

Overall heat exchanger conductance 

Sped f i c vol ume 

Mass flow rate 
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Mass fraction liquid 
Volume fraction liquid 

Z Zeta-factor , flow rastriction factor for nylon fin seals 

G reek 

S Fi n thi ckness 

®HX Heat exchanger effectiveness 

Fin surface efficiency 

TIq Overall finned-tube surface efficiency 

\i Dynamic vi scosi ty 
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MAIN 



Figure A-1 . JTX Computer Program Subroutine Structure 


S-11013 













MAIN PROGRAM 


Call 


HEREADi 


Reads "He 1 i um/Data" H, V, |i, k vs P, T; 

Program internal units are PSIA, °R, BTU, LB, FT, HR 

Read Geometrical Input (inches), (Watts), Z (ZETA-Factor) 

Read Input Conditions: 


T,, P,. T3 

Make Initial Guesses 


"^2 = ■^2MIN = (4.0 °R) 


W = 2.722 X 10 ‘' " 
where AP = P 



Vel 


1/2 


= LB/HR 


f = 0.025 

X 


V = (V2 + \)/ 2 , V2 = f(P^ , T^) 

= f(H2, SAT. L IQ.) 


Call 


NEWTON 


Which calls JTX numerous times, gives all cycle points, W, e 
Calculate: 


HX 


Xl = 

\ - ViLllW “ 3 V 

= W(H^ - h^)/3.4i443 

/ 

*^NET " *^T0T " \eAK 
^JTX ~ ^NET'^^TOT 


Wr? te Output 
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SUBROUTINE JTX; CYCLE ANALYSIS 


Call 

DPIHX 



4P, = f (P,. T,, 
P 2 = P, - AP, 

^ 2 ’ W) 

Call 

EXPDR 



4P^ ’ f (P 2 . Tj, 

W) 

Call 

SATPRP 
P 3 = t (T 3 ) 


Cal 1 

DPOHX 



APq = f (P3, T^, T^, W) 

T3 Is approximated by since is stiH unknown 

= P, - AP 

5 3 0 

First convergence closure 


CLX1 = 


LOG^ (AP^) - LOG^ (Pj - P3>yL0G^ (P^ 
Properties Lookups 


SATPRP 


THERMO 


H4= (P 3 ). (H4= H3^) 

Hi = f (p^, t^) 

»2 = f (P 2 . T^) 


HX Heat Balance 



H 5 = ^ Hi - 

Cal 1 ItHERMO 

^5 “ ^ (•’5’ H3) , I^T^^Ti, Set Flag s = i . O-for printout 

later i n iMAI N( 
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SUBROUTINE JTX (Continued) 
HX Heat Capacities 


CpH = (H, - H2)/(T, - T^) 

Cp|_ = (Hj - H^)/(Tj - T^) 

^MIN ~ ^^PH’ ’'PL^ ’ ^HAX ~ ^‘'PH’ “"PL^ ’ ^ ” 


Cal 1 


HEATX 


UA = f (P,. T,, -Pj. Tj, P3, T3, Tj, W. Cp„, Cp|_) 
HX Effectiveness 

"tu = ^hin) 

*HX = I' - EXP [-\u - C-EXP[-N(„ 0 

if 0 .999999, 


^5 ^4 ■*' ®HX " V 

^5 ^ ^PL " Sin 


PL 


'HX 


CpL (Tj - T^^)/ (T^ - 


Second Convergence Closure 
CLX2 = (T^' - T^)/T^ 


SUBROUTINE NEWTON; CONVERGENCE CRI TER I ON 


Guess T^ and W i n 


MAIN 


Call 


JTX 


To obtain CLX1 and CLX2 


Calculate Error Vector 

ERR = \! (CLXI)^ + (CLX2)^ 

Test Convergence 

If ERR ^ CONV, then T^ and W are okay 
CONV = 0.005 ( 0 . 5 %) 


r /r 
^MIN MA> 
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SUBROUTINE DPIHX 


i_HI GH-PRESSURE-SI DE (INSIDE TUBEl PRESSURE DROP 


Calculate Geometry Parameters (All Dimensions Were inches) 

^NO-FlN ~ (Built in — accounts for inlet section of 

tubing upstream of HX.) 


“h - 



G = W/A, 
Fluid Properties 


TRANSN 


M-1 = f (Pi» ip 

P'2 ~ (^2’ ^2^ ’ *^2 ' ^ approximated by P since P is still 

unknown ^ 

PT = 


THERMO 


Vi = f (P,, T^) 
\ ' f (Pj. V 
w - (V, + V )/2 


Pressure Drop 

Re = Dj^ G/jjf 


f. = 0.029 (Dj/D + 0.304 (Re)'°*^5 


APi = 


^1^4 4q-fin^ 


D. 



— 2 
V G 


288 g 
L =>0 J 


= PSID 


subroutine DPOHX ; LO WrPRESSURE-S I DE (OUTSIDE FINNED TUBE) PRESSURE DROP 
Calculate Geometrical Parameters (All Dimensions Were Inches) 


^ ^ ^ 4 4 4 / 0 ^^ 4 ) 


144 
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p' = P3 - AP„ 

5 o 

^5 ''SMIN’ = '’3 " %iN’ ES] • The output 

value of P^ will then be ”0.10 TORR”, signalling 
that aPq is too large, but program continues. 
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^5 %IN’ jTHERMOl = f {P^ , T^) , v' = (V^^ + v' )/2 

ap' = dpv * 

0 

^ = fj - ''^o 

At this point an Iteration is performed on v/ until P^ and P_ are 

t> 5 5 

converged (0.5^), or 20 iterations are made, then last AP' becomes 

o 

final answer. 

I f < *^5M1 N this iteration, then return to [jTX | as noted 

above, and P^ will be "0.10 TORR". 

SUBROUTINE EXPDR : EXPANDER PRESSURE DROP 

Tv«fo Capillary Tubes, each Divided into 20 Increments of Length, each 
Tube Increment Uses Same Equation for aP 


Re = D .G/fl where i = 1 or 2 

0 f 

0.5 “0.25 

f = 0.029 (D„,/D ) + 0.304 (Re) 


Sudden Expansion Loss Between Tubes 

r . 2 






Exit Expansion Loss 

K^X = 1*0 APg = 288 g 

^c 

Fluid Properties Depend on increment Inlet Pressure, P , and Temperature, T 

X X 3 

Lookups Depend on Pressure Level 

IthermoI 


^2 " ^ ^'^ 2 ’ "^ 2 ^ 
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SUBROUTINE EXPDR (Continued) 


SATPRP 


PgL = f (H2» sat. LIQ.) 

Pj^ = 0.730 PSiA (BUILT IN) 


If P > Pc, , 
X SL 


V = f (P^, H^) 

■'x = ' ('x’ 
tr=f (P^, T^) 


THERMO 


TRANSP 


If P^ < P^ ^ Pgl_, V^, IJ-L = f (P^» SAT. LIQ.) 

Hy, V^, = f (P^, SAT. VAP.) 


X = (H^ - Hl)/(H^ - H^) 


V = + X (V^ - Vj.) 

il = |j,j^ + X (p,^ - lA^) 


Mixed-mean, or Homogenized 
Properti es 


If P^, V, il= f (P^, SAT. VAP.) [SATPRP 

Normal Calculation 

P, - P^P2 - P, - iP, . P3 = Pj - iPj, ETC ... . P„ = P^_, - /iP^., 

Exit of section, P , . = P - aP ; 

n+1 n n 

Exit of expander, P = P , - aP,. - 

’ex n+1 = 1.0 

AP = P^ - P 
X 2 ex 

Conditional Calculation 

Does not stop program, just assumes W is wrong, NEWTON should 
determine next step. When flow is too high, the exit of any 
step could gi.ve negative pressure. Calculate an extrapolated 
(to end of line) exit pressure 

•P'x = P| “ Lj)/(AL)p where L. = i(AL) 


^^x = ^2 
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SUBROUTINE HEATX ; HEAT EXCHANGER OVERALL CONDUCTANCE 


Geometrical Parameters same as [Fp! HX [ and [pPQHX | , with addition of: 


A. = TT D. L^/144 

Afs = (n L^/2p) - D^^)/144 


Fluid Properties same as jPPlHX [ and |pP0HX j , with addition of: 


Pr = Cp p/ k = f (P^, T^, P^, T 2 , P 3 , T 3 , T 3 ) 

Heat Transfer Coefficients ^ ^ 

Inside tube: h. = 0.023 (Re.)°'^^ ^ ^p^^ 

Outside (Fins): h = 0.221 (Re C„, G (Pr 

o 'o' PL o ' 

Fin Efficiency (Applies to Copper Only) 

Call 


KFIN 


^f' ~ ^ ( 7 ), T - (T^ ^ ^2 ^ ^3 ^ 7q)/4 Subscripts: 

f = fin (OFHC CU) 
t = tube (STAINLESS) 


i24 h 


ML^ = 


17 ? 


ML 

■n = TANH f 
’f ML^ 

Overall Conductance (BTU/Hr-°R) 
UA = 


1 LOG (P /P.) 
I ^ e ^ o 

hj A. 2 tt L^ 
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• APPENDIX B 
THERMOPHYSICAL PROPERTIES 



APPENDIX B 


THERMOPHYSICAL PROPERTIES 

DEVELOPMENT OF PRESSURE-ENTHALPY DIAGRAM 
Introduction 


In order to analyze the performance of a Jou I e-Thomson expander operating 
into the he I i urn II (superfluid) region, that Is, at temperatures below 2.177 K 
(saturation pressure of 0.0497 atm), it Is most convenient to utilize a pres- 
sure-versus-enthalpy (P-H) diagram. Since a diagram of this type could not 
be found *n the literature, one was generated from a recent publication of 
hei mm I data (Reference 5), together with several other sources of certain 
helium II data (References 6, 7, and 8), and some key assumptions. The follow- 

I ng paragraphs outline the procedure used to generate the helium P-H Oiaqram 
of F igure B-1 , ^ 

He! ium I Above 1 .0 Psia 

Referring to Figure B-1, all the data above 4.0 R and 1.0 psia are taken 
directly from McCarty, NBS-TN-622 (Reference 5). 

Helium II Coexisting with Vapor 

The region below the lower lambda point, "LLP" in Figure B-1, and to the 
left of the saturated He I vapor line is constructed as follows. The phase 
(P-T) relationship, plotted in Figure B-2, is taken from Cook, p. 324 
(R6f©r6nc9 6), which Is based on +he NBS 1958 He-4 fempera+ure scale* A 
preliminary estimate, of the liquid He II line was calculated from specific 
heat (Cg)* data (Cook, p. 344) and the enthalpy equation; 


H 


L 


T 

/ Cg dT + (P/p) 


(B-1 ) 


where the second term is assumed to be negligible, and Cg « 0 at To = 1.5 R. 

This probably slightly underestimates Cy (heat capacity at constant volume) 
for LHe II. Then a preliminary He I vapor line was obtained by adding the 
latent heat of vaporization (from Cook, p. 348) to Hl* Extrapolating latent 
heat values from McCarty to 3.919 R (lower lambda point as given by McCarty) 
showed a 2.3 percent discrepancy with the earlier published data (Cook, p. 348 
and Gorter, Reference 7, p. 459; same original source). ^However, the enthalpy of 


is the heat capacity measured at the vapor pressure. 
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VAPOR PRESSURE, PS I A 
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Figure B-2. Vapor Pressure of Helium-4 including Helium II Region 
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saturated vapor at 3.919 R was only 1.2 percent different. Therefore, the 
saturated vapor values below the lambda point were adjusted to make a smooth 
curve with McCarty values above 1.0 psfa. Using the latent heat data from 
Cook and Gorter with this saturated vapor line resulted in a shift of the H 
curve. The LHe II enthalpy curve immediately below the lambda point is not 
exactly known (Cook, p. 347). The curve as shown In Figure B-1 near the lamda 
point largely reflects the difference between the extrapolated McCarty and Cook 
latent heat of vaporization data, and the sharp peak in the specific heat curve 
(used in Equation B-1). The liquid He II curve is felt to be reasonable accu- 
rate below 3.8 R (approximately 0.6 psla). 

Low-Pressure Superheated Vapor 

Isotherms below 1.0 psia on the vapor side of the saturation dome are 
constructed from McCarty’s heat capacity data, Cp, extrapolated to the zero 
pressure (Ideal gas) value of 1.241 Btu/lb-R. Then the enthalpy is calculated 
from 


T 

^ = '^SV 


where is the previously calculated saturated He I vapor enthalpy and Tg 
is the saturation temperature. Departure from ideal gas behavior is insigni- 
ficant below about 0.1 psla, hence the isotherms become vertical lines. 

TABULATION OF PROPERTIES 


The following pages present a complete tabulation of the thermophysical 
properties of helium used in the JTX computer program, as well as the thermal 
conductivity for stainless steel and OFHC copper in the temperature range 
below 50 K. 
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